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ABSTRACT 

We present the results of a s urvey for trans-neptun ian objects (TNOs) based on Subaru archival 
images, originally collected by ISheppard et al.1 (|2005[ ) as part of a search for irregular satellites of 
Uranus. The survey region covers 2.8 deg^, centered on Uranus and observed near opposition on two 
adjacent nights. Our survey reaches half its maximum detection efficiency at i?=25.69 ± 0.01. The 
objects detected correspond to 82 TNOs, five Centaurs, and five irregular satellites. We model the 
cumulative number of TNOs brighter than a given apparent magnitude with both a single and double 
power law. The best fit single power law, with one object per square degree at magnitude i?o=22.6_Q'4 

and a slope of Q;=0.51_Q'g, is inconsistent with the results of similar searches with shallower limiting 
magnitudes. The best fit double power law, with a bright-end slope ai=0.7lg'^, a faint-end slope 
a2=0.3to2j a differential number density at i? = 23 cr23=2.0lgg and a magnitude break in the slope 
at i?eij=24.3^Q'^, is more likely than the single power law by a Bayes factor of ~26. This is the first 
survey with sufficient depth and areal coverage to identify the magnitude at which the break occurs 
without relying on the results of other surveys. 

We estimate barycentric distances for the 73 objects that have 24 hr arcs; only two have heliocentric 
distances as large as ^^50 AU. We combine the distribution of observed distances with the size 
distribution that corresponds to a double power law luminosity function to set a tight constraint on 
the existence of a distant TNO population. We can exclude such a population at 60 AU, with 95% 
confidence, assuming it has the same size distribution and albedo as the observed TNOs, if it exceeds 
8% of mass of the observed TNOs. 
Subject headings: Kuiper Belt - Outer Solar System - Trans-neptunian Object 
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1. INTRODUCTION 

The remnants of the protoplanetary disk, now in the 
form of trans-neptunian objects (TNOs), offer a unique 
way to study the evolution of the solar system. The 
TNO size distribution is defined by its initial properties, 
coUisional history, and the formation and evolution of 
the giant planets (iKenvon fc BromlevI 120041 : iPan fc Saril 
I2005t IKenvon et al.l 120071 ). The orbital dynamics of the 
TNOs is largely governed by interactions with Neptune, 
and the radial distribution of T NOs also depends on 
the giant planets' evolution (see iMorbidelli et al.l |2007| 
for a review). It has been suggested that the ra- 
dial extent of TNOs was truncated by a close passage 
of a star d uring the early stages of t he Solar System 
formation (iBrunini fc Fernandez ' 1996; Ida e ^l. 1 120001 : 
iKobavashi fc Idall2001l : IKenvon fc Bromlev.20"oir 

A number of large-scale investigations that will sig- 
nificantly advance our understanding of the outer so- 
lar system are currently being designe d, tested, and ex- 
ecuted. Pan-STARRS (|Jewitti |2003[ ). given its cover- 
age of the sky and time baseline, promises an accu- 
rate determination of the statis tical properties of the 
trans-neptunian regio n. LSST (iTvson fc Ange] 120011) 
and SkyMapper (Kelle r^t al.ll2007l ) will extend the sur- 
veyed sky to the southern hemisphere. The New Hori- 
zons (NH) mission will give unprecedented views of the 
trans-neptunian space by approaching Pluto and other 

^ Based on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan. 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138; cfuentes@cfa.harvard.edu 



TNOs in mid-2015. Nevertheless, there remain impor- 
tant questions that will not be answered by these stud- 
ies. These large synoptic surveys will necessarily have a 
shallow limiting magnitude. Deep surveys like this will 
continue to be the only window into the smallest and far- 
thest objects in the Solar System. The answers to these 
questions can infiuence how these projects are carried 
out and how their resulting data are interpreted. The 
distribution of faint objects will matter when large sur- 
veys choose fields to be covered more deeply. The TNO 
size distribution and radial extent of TNOs are among 
the questions that will rely on pencil-beam surveys to be 
answered. 

Since the discovery of 1992 QBi (jJewitt et al.l Il992l ) 
a number of wide-field surveys for TNOs h ave been 
completed (iJewitt et al.l [199 8: Chi ang fc Brow n 199^ 



Larsen et a l. 2001; ' Truiillo et al.ll200ll:lTnililio fc Brow: 



mg_ 



2noH iMillis et all [20021: iTruiillo fc Brownl 12001 
Elliot et all 120051 : iLarsen et all I2007D . LTaddition 



to determining much of the dynamical structure of the 
trans-neptunian region and identifying large, bright 
TNOs that are amenable to follow up observations, 
these surveys constrain the bright end (i?<24) of the 
cumulative luminosity function of TNOs, the number of 
objects per square degree brighter than a given magni- 
tude. This quantity has consistently been measured to 
be a power law of the form 



S(i?) = io"(^- 



-Ro) 



(1) 



where i?o~23 is the magnitude at which one expects 1 
object per square degree and a is the slope of the distri- 
bution. 
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Ground based efforts have also focused on detecting 
fainter TNOs with deeper imaging of narrow areas of the 
sky. Many ha ve been successfuUy conducted in recent 
vears (, Gladma n_ et al.lfr9 98. 2001; AU en_et al.ll2001l . l2002l : 
iPetit et al, 2004, 2006; Frascr ct al. 2008). These surveys 
have been concentrated near the echptic plane and reach 
limiting magnitudes as faint as R^26. These surveys 
also find that the cumulative surface density of TNOs is 
consistent with a single power law. 

However, the deepest search to date, using the Hub- 
ble Space Telescope with the Advanced Camera for Sur- 
veys and reaching a 50% detection efficiency at mag- 
nitude i?=28.5, found 25 times fewer objects than ex- 
pecte d from extrapolating the brighter (i?<25) distribu- 
tion fBern stein et al.l|200 4). Their result indicates there 
is a break in the cumulative surface density of objects 
n ear R^25. 

iBernstein et al.l ()2004[ ) necessarily relied upon the re- 
sults of other surveys to assess the deviation of the cumu- 
lative density of objects from a single power law over a 
range of magnitudes. However, it is difficult to combine 
the results of different surveys to obtain a well-calibrated 
sample of the trans-neptunian population. 

Dynamical biases in latitude and longitude, change the 
local density of objects and the relative abundances of 
excited and classical objects depending on the direction 
in which a survey is conducted. This can be seen in 
the variety of results found in the literature; a nice re- 
analysis and summ ary of some surveys is presented by 
iFraser et al.l 120081 . Different surveys sample a variety 
of ecliptic latitudes and longitudes, use various analysis 
methods, or vary in observing conditions. It is necessary 
to determine and correct for the effects of these differ- 
ences to characterize the physical and dynamical prop- 
erties of the TNO population. For bright TNOs, large 
synoptic surveys will determine many of the biases in the 
observations as well as in the population itself. However, 
for fainter TNOs, the simplest way to overcome these 
difficulties is to observe a single region of the sky. 

All these surveys use the "digital tracking" method 
JGladman & Kavelaars 1997; Gladman ct al. 1998, 2001; 
[Alien et al.l 1200 iL |2002) , by which a series of consecu- 
tive short exposures are digitally shifted and coadded to 
match the apparent motion of real objects. That method 
has proven very useful in improving the sensitivity of 
these ground and space based observations. However, 
this method relies on how fine the grid of velocities sam- 
pled is, the good quality of a template image to subtract 
from each exposure, extra processing of the images and 
a trained operator to filter false positive detections due 
to saturated stars or other artifacts. Our results were 
obtained by linking detections in three different images 
described in SjS] Our method's data reduction is more 
direct, requires less human interaction and is easier to 
photometrically calibrate. 

The radial extent of the classical TNO population is 
not known with certainty. Although there is evidence for 
a sudden decrease in density at r '^ 50 AU |Truiillo ct al. 
I2OOII : iGladman et al.l [2001h . the existence of a second 
farther population near the ecliptic is difficult to rule out, 
due to the bias against detection of more distant, fainter 
objects. We are slightly more sensible to distant, slower 
moving objects. Since we do not rely on the construction 
of a template field, usually made with data taken on the 



same night, that increases the noise and would subtract 
signal from very slow movers. 

The objectives of this work are to better constrain the 
expected break in the TNO luminosity function using a 
single survey and to better understand the lack of detec- 
tions at large heliocentric distances. In the next section 
we describe the data and the processing of images. In §[3] 
we present our moving object detection algorithm. We 
discuss the control population and detection efficiency of 
our method in § [H In the final two section we present 
the results of our survey and disucss their implications 
for the size and distance distribution of the TNO popu- 
lation. 

2. DATA 

The observations considered in this project were 
taken on UT 2003 Augus t 29 and 30 with Suprime- 
Cam ("Mivazaki et al.|[2002f ) mounted on the Subaru tele- 
scope. Suprime-Cam is a mosaic camera with 10 CCDs, 
each with 2048 x 4096 pixels. Each mosaic image has a 
field of view of 34' x 27'. We used SMOKA, the electronic 
archive of the Subaru Telescope (jlchi kawa 2002), to re- 
trieve observations taken in August 2003 in the vicinity 
of Uranus, near opposition. The fields were originally 
observed bv lSheppard et all ()2005( ). They surveyed a to- 
tal of 14 fields, with an areal coverage of 3.57 deg^ over 
the course of two nights. All exposures were taken with 
the "Cousins R" red filter, well-matched to the colors of 
outer solar system objects. 

The objective of the original inves tigation was to 
discover uranian irregular satellites. ISheppard et al.l 
(2003) recovered all previously known uranian irregular 
satellites and discovered two new such satellites. The 
faintest satellites dete c ted have magnitudes at R~25.5 
(jSheppard et all 120051: [Kavelaars et all I2004D . On the 
first night, the observers took two or three exposures of 
^7 min of each field, separated by half an hour on the 
first night. They re-observed those fields with two expo- 
sures during the second night, with the pointings shifted 
to maintain the same positions relative to Uranus. The 
survey was designed to discover objects during the first 
night and to obtain better orbital information using the 
second night's data. 

We chose this particular data set for the following rea- 
sons. The data set is sensitive to B < 25.5 mag nitude 
object s, in the magnitude range in which Bernstein et al.l 
(|2004f ) find the break in the TNO cumulative function to 
be. This sensitivity is reached in a single exposure, avoid- 
ing the difficulties associated with combining different 
images. There are 11 fields (2.8 deg^) with 3 exposures 
on the first night, permitting a simple search for moving 
objects. The fields were observed very close to opposi- 
tion, allowing a reliable distance estimate from the rate of 
motion with only a 24 hr arc. The sky coverage is large 
enough to expect the discovery of ~100 TNOs, allow- 
ing a significant constraint on the cumulative luminosity 
function. Finally, the data were easily obtained from the 
SMOKA system, after the 18 month proprietary period. 

We performed the usual calibration of the images. For 
every image, we performed an overscan correction, trim- 
ming, bias frame subtraction, and fiat-field division using 
standard IRAF"^ routines. Calibration frames taken dur- 

^ IRAF is distributed by the National Optical Astronomy Ob- 
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ing these observations were obtained from SMOKA. 

3. MOVING OBJECT DETECTION 

The apparent motion of outer solar system objects 
viewed near opposition is primarily due to the Earth's 
translation. For objects at the distance of Uranus the 
apparent motion can be as large as ~6 "/hr. For TNOs 
this rate is typically ~3 "/hr. This motion, with respect 
to background stars, is readily detected even in the short 
(~1 hour) time baseline of this dataset. 

To find TNOs, Centaurs, and irregular satellites 
in this data set, we use a variant of the search 
algori thm described and implemented by iPetit et ahl 
(|2004f ). This method is simi lar to that used in other 
TNO surveys (for example. iLevison fc DuncanI 199CIJ : 



Irwin et al . 1995: J ewitt fc Luulll9 95: Truiillo ct al. 2001'; 
Millis et al...2002. ). The algorithm detects moving objects 
by comparing the positions of all point sources in each of 
three images of a patch of sky taken in the same night. 
Thus, as mentioned earlier, we restricted our search to 
the 11 fields for which there were 3 images taken on the 
first night. The individual steps in the algorithm are as 
follows. 

First, for each search field we determine an astrometric 
solution for the first image of the night. These astromet- 
ric solutions are used later to guide the insertion of syn- 
thetic m oving objects. We used the 2MASS point source 
catalog (jCutri et al.1 l2003f ) as an astrometric reference. 
The RMS in the astrometric solution was typically of 
0.2 " or lower (close to the catalog's precision). The rel- 
ative errors on the astrometric solutions for both nights 
were comparable to the tolerance of the search algorithm. 

We then register the second and third images with the 
first image of the night. This allows for very accurate 
positioning of stellar-like objects with respect to each 
other. This is done for the individual CCDs, rather than 
for the entire mosaic. The successive CCDs images are 
linearly interpolated, automatically, to the first using the 
positions of the back ground stars and routines available 
in the ISIS package () Alardl [2000l ) . When thse routines 
failed to converge (due to numerous bad pixels or satu- 
rated stars) , we align the images interactively using rou- 
tines from IRAF. 

At this stage, we insert the population of synthetic 
objects that will be used to determine the detection effi- 
ciency of the search, as described in § |4l 

We then use two different algorithms to search for point 
sources. The first of these is a wavelet transform source 
detection routine (see lPetit et al.|[200^ for a description). 
The second is the pubhcly available SExtractor pack- 
age (|Bertin fc Arnoutill996f ). which calcuates the local 
image background RMS, convolves the image with a user- 
specified kernel, and then identifies groups of pixels with 
values exceeding the background variation by a given 
value. These two approaches have very different false 
detection characteristics. Thus, we consider the inter- 
section of detections from both routines. (We use the 
flux information given by SExtractor for the photome- 
try described in § [5]). We use a detection threshold of 
2.6ct. For SExtractor this corresponds to four adjacent 

servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



pixels with values that are at least 1.3 times of the lo- 
cal background variation. At the end of this stage, there 
are three lists of sources, one list for each image. This 
results in up to ~50,000 detection in each mosaic im- 
age. Note that the expected motion of trans-neptunian 
objects during a single exposure ('~7 minutes) is small 
compared to the typical FWHM (0.7 "). Thus, trailing 
does not significantly affect the source detection. 

In order to identify moving objects among all the point 
source detected, we apply a series of filters that eliminate 
individual detections, as well as sequences of detections, 
that are not consistent with the TNO population. 

We first reject all detections that corresponded to sta- 
tionary objects, i.e. stars and galaxies. For each list of 
detections, we eliminate those for which there is a cor- 
responding detection within 0.05 " in at least one of the 
other two lists. We deliberately chose a small threshold 
in order to not diminish our sensitivity to very slow mov- 
ing TNOs. This stage typically reduces the number of 
detections to '^10,000 per field. 

The next step is to search for linear motion among the 
non-stationary detections. We identify all groups of three 
detections in the successive images that are consistent 
with straight line motione (within 15° of the ecliptic), 
with a constant angular rate between 0.5 and 10 "/hr. 
The parameter space is chosen to include the expected 
rate and direction of TNOs. We consider all combina- 
tions of detections in the three different exposures whose 
fit to a line had an RMS of 0.3 " or less. These criteria 
are met by ~1,600 combinations per field, nearly all of 
which are synthetic TNOs (see 35). 

In the final stage, the search program outputs an im- 
age with all the combinations found, showing a stamp 
centered on every detection. We visually inspect these 
images to accept or reject a given object, images. This 
method allows the spurious and acceptable detections to 
be rapidly distinguished ('^30 min/field). Typically ~20 
objects are rejected in this stage per field, the majority 
being optical artifacts, bad pixels, extended objects or 
some combination of the above. 

4. CONTROL POPULATION AND DETECTION 
EFFICIENCY 

Since our observations are flux-limited it necessary to 
account for detection biases when estimating the intrin- 
sic number of TNOs as a function of magnitude. We 
characterize our search using a population of synthetic 
TNOs inserted just after the images have been calibrated 
and their astrometric solutions determined. The pro- 
cedure is done for each mosaic field, rather than CCD 
by CCD. This process nicely accounts for the possible 
motion across detector boundaries. The same synthetic 
populations were used for the s econd night. 

We used the Orbfit routines (jBernstein fc Khushalanil 
l2000f l to create a realistic population of synthetic TNOs. 
The characteristics of the population were chosen to span 
the range of observational properties expected of the 
TNOs. The position of an object on the sky at the time of 
the first exposure was drawn from a uniform distribution 
that encompassed the FOV of the mosaic. Objects were 
implanted with distances between 20 AU and 200 AU, 
or alternatively 0.7 "/hr to 5.0 "/hr. The proper motion 
and radial velocity given to the object are taken from 
a distribution that encompasses the possible rates found 
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in the Solar System. This initial position and velocity 
vectors are only accepted if they correspond to a bound 
orbit. If so, we use the Orbfit routines to calculate the 
RA and Dec of the object at the beginning and end times 
of each exposure. We translate these sky positions into 
locations on the mosaic using the astrometric solution 
derived earlier. 

For each exposure we compute a model for the PSF. 
The model is the average of '^10 bright, isolated stars for 
every CCD. Given the known magnitude of the synthetic 
TNO and the measured zero point, and accounting for 
transparency changes through changes in the flux mea- 
sured in the PSF stars, we use IRAF routine to insert 
PSFs with this flux at the calculated positions. We in- 
serted objects from 22.5 to 26.5 mag, which spans the 
magnitudes of the TNO population we expect to find. 
The flux of each object includes photon noise. We did 
not consider variable objects, as this is unlikely to be 
significant on ^^1 hr time scales. We include the effect 
of trailing by dividing the flux among several PSFs in- 
serted at positions linearly interpolated between those 
at the beginning and end of the exposures. This process 
takes into account any background, transparency, see- 
ing, and focus variations that might affect the limiting 
magnitude. Using this PSF model from each image, we 
implant a set of ~2,000 objects per fleld. This results 
in a sufficient number of synthetic objects per CCD to 
sample the detection efficiency as a function of position. 

Since we are counting objects up to a certain brightness 
and our model describes the underlying TNO population, 
it is essential to estimate what fraction of the population 
we detect as a function of magnitude. In Figure [T] we 
include a histogram of the fraction of objects that were 
recovered in each magnitude bin. We implanted 25,074 
objects in 11 fields, recovering 17,195 of them. 

When plotting the cumulative function we used the lo- 
cal efficiency function, each detection is weighed by the 
number of objects recovered in the same field and within 
the observational magnitude error. The detection effi- 
ciency could vary from field to field. Since all fields were 
taken in the vicinity of Uranus, efficiency could depend 
on location. However, its statistical effect was negligible 
on the efficiency. 

For the statistical analysis the effective efficiency func- 
tion will need to be integrated. Since it is simpler to inte- 
grate analytical expressions, we used the tot al efficiency 
functi on, that considers all fields. Following iPetit et al.l 
()2006( l. we represent it by 



viR) = J I 1 



tanh 



R — i?5o 



1 — tanh 



R-R 



50 



(2) 
where the best fit values are A=0.88±0.01, i?5o=25.69± 
0.01, u;i=0.28 ± 0.04 and w;2=0.88 ± 0.15. The errors 
were obtained with a Markov Chain Monte Carlo simu- 
lation. The parameter A corresponds to the maximum 
efficiency, achieved for bright objects. R^q corresponds 
to the magnitude at which the detection efficiency drops 
to half the maximum values. The parameters wi and 
W2 characterize the abruptness of the decline of the de- 
tection efficiency. Figure [1] shows the average efficiency 
function for our data set. 

The efficiency could also depend on the rate of motion. 
We construct a rate-analog to the magnitude efficiency 



(See Fig. [21). The detection efficiency is nearly indepe- 
dent of rate, but our method is slightly less efficient at 
larger rates. A faster moving object that is detected in 
the first image has a greater chance of falling close to 
a background star, or moving outside the field of view, 
thus the detection efficiency declines with the rate of mo- 
tion. The lowest bin plot in Fig. [2] is 1.5 "/hr. Since 
we implanted objects to have a population with a con- 
stant surface density that bin is not well sampled. Even 
though we were able to recover objects planted with rates 
as slow as 0.7 "/hr (parallax for objects at 200 AU) we 
consider a more conservative limit. The rate at which an 
object moves 1-FWHM in 45 minutes, the shortest sep- 
aration between the first and third exposure, is 0.9 "/hr 
(150 AU). 

To properly account for detection biases, both real and 
control objects must go through exactly the same vali- 
dation procedure. We did not unveil the fake object list 
until all objects were recognized as moving objects, either 
real or planted. 

5. RESULTS AND ANALYSIS 

We found 92 moving objects, five corresponding 
to known irregular u rani an satellites (those f ound by 
iKavelaars eFal] 120041 and lSheppard et al.l [20051) . five to 
Centaurs, and 82 to TNOs. The satellites that were 
missed were blended with stars in one of the images and 
hence were not found by our algorithm. 

We present our detections in Table [TJ For each TNO, 
we list its internal designation, its position at the time of 
the first exposure (also listed), and its estimated magni- 
tude with uncertainties (along with an independent es- 
timate of the photometric uncertainty). We also list the 
measured sky plane rates of motion of the TNO, two es- 
timates of the distance to the TNO (one suited for two 
night's data and another based only on parallax, both 
explained later), its orbital inclination, and separation 
from Uranus at the tim e of disco very. 

Three standard stars ([Landolti 19921 were used to ob- 
tain the zero point and airmass dependence of the pho- 
tometry. These were PG2213-006C (V=15.11 ± 0.0045, 
V-R=0.426 ± 0.0023) SA-92-417 (V=15.92 ± 0.0127, V- 
R=0.351 ± 0.0151) and SA-92-347 (V=15.75 ± 0.0255, V- 
R=0.339±0.0295 ). Since their colors a re similar to those 
of typical TNOs ([Peixinho et al.ll2004D we did not apply 
a color correction. We checked both nights were photo- 
metric and stable. The possible dependence on seeing 
(FWHM) was also investigated, finding it to be unim- 
portant. The correction term was negligible compared 
to the airmass correction. Every detection's magnitude 
is calculated, using the following formula: 

i? = 27.36 -2.5 log /s/i- 0.09 a:, (3) 

where /s corresponds to the flux in a 5-pixel aperture, t is 
the time in seconds and X is the airmass. This equation 
accounts for an average 0.34 mag aperture correction be- 
tween the known magnitude of a synthetic object and its 
magnitude measured with a 5-pixel aperture. The search 
algorithm requires an object to be found in all three expo- 
sures giving three independent magnitude measurements 
that we average to obtain the results shown in Table [TJ 
The errors given on the magnitude values correspond to 
the error on the flux. 

In Figure [31 we plot the photometric errors, showing 
them to be '^0.1 mag. The magnitude dependence of 
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the uncertainty is shown in figure [4l We estimate the 
uncertainties empirically, calculating the standard devi- 
ation as a function of magnitude. We then fit a second 
degree polynomial, overlayed in Figure ID This estimate 
is shown for each real object as ^Rmag in Table [TJ 

We can accurately approximate the apparent motion 
of a TNO over 24 hours as a straight line with a con- 
stant rate. We include the measured right ascension and 
declination rates in Table [1] The apparent motion of our 
objects compared with that which parallels the ecliptic 
is plotted in Figure [5] 

Near opposition, the change in the rate or direction 
of motion over 24 hours is negligible, making it easy 
to predict where the real objects would be on the sec- 
ond night. However, nine objects were not found on the 
second night. Our method is only ~ 90% efficient for 
the brightest objects on the first night, with 10% lost to 
blending with field stars. It is expected that more than 
10% of the objects will be lost on the second night, be- 
cause of confusion with stationary sources and because 
they are more likely to move outside the field of view 
over 24 hours. 

We used the observations on the second night to im- 
pro ye the distance determination w hen possible. We use 
the iBernstein fc Kh ushalanil ()2000D Orb fit routines to es- 
timate plausible orbital elements assuming there's no ac- 
celeration in the direction tangential to the plane of the 
sky. For a 24-hout arc, this results in a '^7% accuracy in 
the barycentric distance {dtari)- For a single night ob- 
servation of objects the error on the distance could be 
unbound. However, since the observations were taken 
near opposition we are able to readily estimate heliocen- 
tric distances (dpar) from the "parallactic motion". We 
assume that the observations are taken exactly at op- 
position and that the orbits are circular. This distance 
estimate is not as reliable as dbari but it serves as a con- 
sistency check. 

5.1. Statistical Analysis 

The probability of our data {D) given a model for 
the intrinsic population (M) is denoted P{D\M)—L{M) , 
where L is the likelihood function. We consider the 
data in our survey as a collection of N detections with 
measu red magnitudes. As derived in ISche chter fc Presa 
(|1976f ) if g{m)dm is the expected number of detections 
between m and m + dm, then the likelihood of a set rrii 
where i=l, • ■ ■ , A'^ is: 



N 



L{M) — exp[— I g{in)dm] I I g{rai)dm. 



(4) 



We are interested i n cha racterizing c/(m). As described 
in IBernstein et all ([2004), we can think of g{m) as being 
the probability of detecting an object and assigning it 
a magnitude m given the survey characteristics and the 
real distribution of objects on the sky. We consider an 
intrinsic differential surface density of objects a that only 
depends on magnitude and is constant over the observed 
area as the model M. For a survey with an efficiency 
function rj, a function of magnitude only, we can write 

oo 

J g{m)dm = Q J 'i]{m)a(m)dm, where 51 is the solid 

— oo 

angle of the survey. 



The likelihood of a model for the differential surface 
density a{m) is then given by: 



_^ — ilf ri{m)a[va)dra 



n / k{m),im)aim)dm. (5) 



L{a) 



This is the probability of finding each object in the set 
of observations at its measured magnitude, scaled by the 
probability of not finding anything else. The function 
li{m) is the probability an object is given a magnitude 
rrii given its intrinsic magnitude is m. 

If we consider the efficiency function and model it as 
relatively linear over the magnitude uncertainty of an 
observation we can approximate our likelihood function 
as follows: 



L{a) 



— p-^'>- S 'n{x)ry(x)dx 



W_ri{rni)a{mi) 



(6) 



This is extremely useful when dealing with a large num- 
ber of objects and surveys. We compared the behavior 
of both exact and approximate likelihood functions with 
our data and found no noticeable differences. 

If we want to sample the likelihood function over its pa- 
rameter space or calculate the total likelihood of a model 
we need to consider priors. That is, the probability of a 
parameter q given a certain model M, P{q\M). These 
priors reflect any knowledge we have over the value of a 
parameter previous to our survey. We chose priors that 
reflect the least previous knowledge into the analysis. We 
chose uniform functions between limits set by our survey, 
indicating our ignorance of those parameters. The total 
probability of a model is: 



C{a)= / P{q\a)L{a,q)dq 



(7) 



We can compare two competing models using their total 
likelihoods by computing the odds ratio: 



021 = 



P{a2\D) P{a2)P{D\a2) £(^2) 



(8) 



P{<Ti\D) P{cji)P{D\cji) C{cji) 

The last equality holds if we do not have a good reason 
to prefer "a priori" any of the two models. The ratio of 
the total likelihoods is called Bayes factor. 

5.2. Single Power Law Model 

One of our goals is to determine whether the results of 
our survey indicate that the cumulative surface density 
can be modeled by a single power law distribution (SPL) 
or if the data favor a more complicated model. We use a 
likelihood analysis to investigate this. 

The likelihood function is related to both the detec- 
tion efficiency of the survey and the differential surface 
density a{R). The important observation for the anal- 
ysis is the number of objects we detect brighter than a 
given magnitude, namely the cumulative surface density: 
E(i?)=/_^ a{x)dx. We use the likelihood function given 

by eq. [5l with 57=2.83 deg^ and r]{R) given by eq. [2l 
For every object we model its photometric uncertainty 

using the analytical model we considered previously, a 

gaussian [li) around its measured magnitude (see figure 

11. 

The single power law model is written as follows: 



ui{R,a,Ro) = aln(10)10"(-"~-"°) 



(9) 
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In Figure [6] we plot the SPL likelihood as a function 
of Ro and a. The previously accept ed values for the 
SPL parameters (a=0.76, i?o=23.3) (jPetit et all 120061) 
are in strong disagreement with our data, lying well out- 
side our 3-cr confidence region. Most of the surveys 
that have consistently measured a slope of q;~0.7 for 
the cu mulative distribution have brighter lim i ting m agni- 
tudes (iGladman et al.lll998HPetit et al.|[200l|2006D. The 



likelihood function (Eq. [5]) with: 



fT2(i?)=C 



^Q~ai{R-23) j^ -[^Q(Q2-Ql)(flc5-23)-Q2(fl-23) 



except ions are IGladman et al.1 (|200in and iFraser et al.1 
(|2008f ). who quote a magnitude hm it of R =25.9 and 
i?=25.6 respectively. iBernstein et al.l ()2004D performed 
a search complete to R=28.5. They discovered far too 
few objects to be consistent with a SPL. 

We check that our bright end sample is consistent with 
the previous surveys with shallower limiting magnitudes. 
In Figure [7] we plot our sample's likelihood function af- 
ter imposing an artificial efficiency limit at i? =24.5. The 
power law index is clearly consistent with the lPetit et al.l 
(|2006( ) result and it shows our sample does not deviate 
from the SPL behavior observed by others for magni- 
tudes brighter than i?~24.5. 

To show that the deviation from a SPL at fainter mag- 
nitudes is not an artifact of our efficiency function, we 
repeated the experiment but instead imposed an artificial 
break at i?=25.2, where our survey is 70% as efficient as 
its maximu m efficiency. The re sult can be seen in Fig. [51 
it shows the iPetit et all (|2006[ ) result is rejected at the 
2-(7 level. 



5.3. Double Power Law Model 

Now that we have shown that our results are not well 
modeled by an SPL, we test a more compficated model. 
Any model that includes a break in the surface density 
distribution will have more free parameters than an SPL. 
Al ternatives with three a nd four parameters were tried 
bv IBernstein et aLl (j2004l ) to explain the aforementioned 
under-abundance of detections. We will focus on the 
"double power law" (DPL) model, the harmonic mean 
of two different power laws. Though a model with three 
parameters would be easier to implement, it does not pro- 
vide the immediate insight into the TNO population that 
the DPL provides. The DPL has four free parameters, al- 
lowing two different asymptotic power law behaviors for 
the distribution (that can be linked to the size distribu- 
tion of small and large objects), a break in the luminosity 
distribution, and a differential density constant. 

The larger number of parameters makes the likelihood 
function more difficult to sample, thus we use a Markov 
Chain Monte Carlo (MCMC) approach (for an MCMC 
review see iTegmark et al.l l2004l ) . We use a Metropolis- 
Hastings algorithm to sample the likelihood function 
with a gaussian proposal distribution. The parameters 
were set to yield a ~25% acceptance rate. We considered 
a run of 100,000 iterations. To check for consistency we 
tried different initial conditions and compared the re- 
sults, no disagreement was found. We also checked the 
performance of our MCMC code with the SPL model. 
In Fig. |S] we show the marginalized probability for both 
parameters a and Rq from MCMC and the exact result. 
There is evident agreement between the two approaches. 

The DPL likelihood function is obtained by replacing 
the corresponding surface number density (Eq. llOp in the 



C-a23(l + 10("^-"i)(-«'=''-23)) (10) 

In Figure [TOj we show the DPL likelihood as a func- 
tion of the bright-end slope ai, the faint-end slope q;2, 
the value of the surface number density at i? = 23 1723 
and the break magnitude Req- All parameters but ai 
are well constrained by the data. Given the small num- 
ber of bright TNOs detected in our survey, the limited 
constraint on ai is not surprising. 

5.4. Cumulative Number Density 

Using the detection efficiency (Eq. [2]) we can estimate 
the number of objects we missed for each object found. 
We construct a cumulative function of the unbiased pop- 
ulation plotting each object individually, representing 
with its detection a number of objects with similar mag- 
nitudes. Since we are plotting a cumulative function, the 
errors are correlated (See Fig. ITT]) . 

We go on to compare the total likelihood of both mod- 
els, as described in § 15.11 A simple way of doing this is 
to examine the goodness-of-fit of the cumulative number 
density. Figure [Tl] shows the data and the best solution 
for the single and double power law cases. Note that 
those power laws correspond to the cumulative number 
densities, Ei(i?)=10"(-«--«°) and ^2iR)^j":^(T2{x)dx. 

It is expected that a DPL gives a better fit to the data 
than a SPL model. The question is whether this bet- 
ter fit overcomes the increased complexity in the model. 
This can be answered calculating the quotient of the to- 
tal bayesian probabilities of the models (Bayes factor, de- 
tails in Appendix iS.ip . If the total probability for a given 
model is larger than another then it is preferred. Using 
the results of the MCMC simulations we compute this 
factor. The resulting total probabilities depend on suit- 
able priors, that reflect our ignorance on the parameters. 
We selected uniform priors for all our variables. For the 
SPL we chose a e [0.35, 0.85], Rq e [21.0, 24.0], while the 
DPL priors were uniform, ai e [0.5,1.0], 02 e [0.1,0.7], 
(723 e [0.5,5.0], Req e [23.0,26.0]. The calculated Oc- 
cam's factor is 0^^=26, meaning that a DPL model is 
more likely to be a better representation for the bright- 
ness distribution of our data. 

5.5. Other Surveys 

IBernstein et al.l ()2004[ ) com bined the results o f thei r 
HST survey with those of Chiang & Brown fl999|}; 
Uladman et al. (1998); Allen et al. (2002); Truiillo ct at] 
((2OOID; iLarsen et al.1 (|200lh : iTruiillo & Brownl (|2003h . 
We include most of the objects listed in that work 
and those conducted since. Table |2| differs from 
IBernstein et al.l ()2004L Table 2) in the exclusion of the 
two w idest searches and the inclusion of tw o newer sur- 
veys (jPetit et al.l 120061 : iFraser et al] 120081 ). as well as 
ours. We excluded the two surveys because of the com- 
plexity in establishing the searche d area near the e c liptic . 
For the sake of comparison with IBernstein et al.l ()2004l ) 
we use the same criteria regarding detected objects as 
well as the caveats provided therein. We included sur- 
veys for which the location of the searched area, effec- 
tive area of the search, magnitude at which the efficiency 
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drops by 50% must be given. We include objects that 
have an observed magnitude where their efficiency func- 
tion is more than 15% the maximum efficiency of the 
survey. We point out that all our detections satisfy this 
requirement. 

We are interested in computing the likelihood of a 
model given the data from each survey. For this we only 
need the list of objects that meet our criteria, an estimate 
of the efficiency function, the surveyed area, and a way 
to translate all measurements to the red filter R for each 
survey. We use the approximation given in (Eq. [6]). Fig- 
ure[T2]shows the 333 objects that we considered. It shows 
the existence of a very pronounced lack of detections at 
faint magnitudes. Our likelihood analysis is summarized 
in Fig. [13] with 1-crconfidence limits for the parameters. 

An interesting aspect of our search is that the data 
has been available since August 2003. Our survey's most 
likely distribution expects '^12 detections for the HST 
field while 3 were found. This provides independent sup- 
port to the existence of a break in the TNO luminosity 
function. 

5.6. Classical & Excited Population 

We use the criteria in lBernstein et all (|2004f ) to identify 
"Classical" and "Excited" objects. TNOs with distance 
at discovery d between 38 AU and 55 AU and inclination 
i < 5 deg are considered "Classical" and the rest are 
considered "Excited" . In Table[2]we list each survey with 
the corresponding number of TNOs in each category. 

This survey was considered by itself and together with 
the surveys in Table [2 We investigated how does the 
DPL luminosity function change when applied to the 
different populations. We repeated the MCMC anal- 
ysis for both populations and for our survey and the 
combined survey. We also considered the priors used in 
iBernstein et al.l ()2004l ). —0.5 < ai,a2 < 1.5 to constrain 
the parameter space. 

The results of the MCMC simulations are summa- 
rized in Table [3l These resu lts are very similar to 
those bv IBernstein et all ()2004[ ). However, we have in- 
cluded three new surveys (This s urvey and those by 
iPetit et al.ll2006HFraser et al.l l2008D. two of which (This 
survey, iFraser et al.l l200 8f) sample magnitudes fainter 
than it!^25.5 where excited objects were specially under- 
sampled. 

5.7. Size, Distance & Inclination Distribution 

The size distribution is closely related to the distribu- 
tion of apparent magnitudes. It is customary to assume 
all objects are located at the same distance and that 
the size distribution is a single power law and hence the 
cumulative brightness distribution is also a power law. 
The parameters of the two distributions are related by 
q = 5a -I- 1, where q is the exponent of the differential 
size distribution (dn — D^'^dD) and a is the exponent 
of the SPL cumulative luminosity function. 

With our rough distance estimates and assuming a 4% 
albedo for TNOs, we can compute the real size distribu- 
tion of the objects in our survey (we adopt rrifj — —27.6 
for the R band magnitude of the Sun). In Figure [14] 
we show the cumulative size distribution for our survey. 
However, the typical error in distance ~7% translates 
into a 0.3 magnitude photometric error, triple the median 



photometric error in our survey (see Figure[3]). Thus, in- 
stead of repeating the statistical analysis for the size dis- 
tribution directly, we transform our luminosity function 
into a size distribution assuming all objects are located 
at 42 AU. The best DPL fits for the luminosity func- 
tion are plotted as a function of size. The solid line is 
the fit to this survey and the dashed line corresponds to 
the fit to the surveys in Table [2] We also consider a toy 
model based on the DPL; it corresponds to two power 
laws with index qi = 5ai -I- 1 and (72 = 5q;2 -I- 1 that are 
joined at the size for which an object at 42 AU would be 
observed to have magnitude Req- We plot the cumula- 
tive function of the toy models for both DPLs to show 
the asymptotic behaviors as a light solid line and a light 
dashed line respectively, both are arbitrarily offset verti- 
cally for clarity. In Figure [14] there is a clear agreement 
between the real size distribution and the fit for the DPL 
models indicates that the assumption that all TNOs are 
at the same distance is justified. 

In Figure [15] we plot the distance and magnitude for 
each object. The distance corresponds to dhari in Table]!] 
with the exception of those objects that were not recov- 
ered on the second night for which we plot the circular 
orbit approximation (dpar)- We consider only the sub- 
set of 73 objects with 24 hour arcs data, with a distance 
error of 5%. 

All but two objects are located at less than 50 AU from 
the Sun, although we are able to detect D = 250 km 
TNOs at distances of 80 AU, with 50% efficiency. This 
lack of distant detections has been noted previously 
(I Allen et all |200H iTruiillo et all 1200 ll : [Bernstein ct al] 
l2004f ) with the recurrent hint that there is an "edge" 
to the Kuiper Belt. 

Given the size distribution that corresponds to our best 
fit luminosity function we are able to calculate the dis- 
tance bias in our sample and obtain the real distance 
distrib ution. We follow the approach of iTrujillo et al.l 
([2001h . The true and observed distributions are related 
by f{r)dr oc P{r)fo{r)dr, where /3{r)^^ ~ J ^ n{D)dD 
is the bias factor and n{D) is the TNO size distribution. 
This is done for 10 magnitude bins between 22nd and 
26th magnitude and independent estimates of the bias 
function are obtained. We used the average to test the 
effect of the DPL size distribution to the distance distri- 
bution of objects, as shown in Fig. [T7] We see an abrupt 
drop in the abundance of objects at r~47 AU, regard- 
less of the size di stribution considered, as has been de - 
scribed by others (jTrujillo et al.ll200lHPetit et al.ll2006[ ). 
However, a DPL size distribution gives a much tighter 
constraint on the existence of a distant population. This 
is due to its much shallower size distribution for small 
bodies as can be seen in the bias correction for the DPL 
for our survey and the one for all surveys. 

Given the fact that we detect no objects farther than 
50 AU we can constrain the surface density E of a dif- 
ferent population located outside 50 AU. At 95% confi- 
dence level, the detection of no objects is consistent with 
an expectation of 3 detected objects. We calculated this 
for the observed population N^xp — ^/n Tj{x)Y,{x)dx, 
where 77 is the detection efficiency of our survey. We 
assume for simplicity that the size distribution of the 
distant population is the same as that we have measured 
for the objects in our survey. We will also assume that 
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each object in the population is shifted to larger heliocen- 
tric distances by the same factor. It is usefiil to define 
the limit on a distant population at distance d as the 
maximum fraction of the observed population's surface 
density that a population can have to be consistent with 
no detections. We denote this fractio n as g{d), following 
the notation in lBernstein et al.l ()2004[ ). For 60 AUwe find 
n — 0. 08, compared to g = 0.17 found bv lBernstein et all 
(|2004f ). Our survey rejects another population with the 
same mass closer than 110 AU. Thus, we place a tight 
limit on the existanc e of a distant populati on. We sup- 
port the conclusion of iBernstein et ahl (j2004f ) that if such 
a population exists, it is either substantially less massive 
than the observed classical Kuiper belt or it is comprised 
of small bodies that are beyond our detection threshold. 
Using the inclination information in Table [T] we can 
show the inclination distribution for the objects in our 
survey. The results are shown in Figu re 1181 This is very 
similar to the results in iBrowia (|2001l ). 

5.8. Mass 

We use the results of our MCMC analysis to estimate 
the total mass of TNOs to which our survey is sensi- 
tive. At each step in the MCMC runs, we compute the 
mass that corresponds to the DPL parameters (again, 
assuming a heliocentric distance of 42 AU and a geomet- 
ric albedo_o£_0;04) . ^^^foUow the parametrization used 
in IBernstein et al.l (|2004l ): 
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where M23 = 6.3 x lO^^kg = 1.055 x lO-^Af® and / is 
the fraction of objects from the given population that are 
located within Q. 

We consider the complete TNO population and the 
Classical and Excited sub-samples. The DPL size dis- 
tribution allows us to compute the value of the integral 
in Equation lll[ however the total mass of a given pop- 
ulation depends heavily on the mean values of the as- 
sumed physical parameters. The mass probability dis- 
tribution is calculated assuming all other parameters are 
fixed. The uncertainties on the rest of the parameters 
(density, albedo, distance and fraction in the surveyed 
area) can be accounted for independently. We consid- 
ered an effective area of ±3 deg from the ecliptic, giving 
ri=21,600 deg^ and that all objects in each population 
are located within that area (/=1). We have also as- 
sumed mean albedo p — 0.04, distance d = 42 AU and 
density p = 1000 kg m^'^. 

In Figure [16] the mass distribution is plotted for our 
survey alone (solid lines) and for the combination of 
all the surveys listed in Table [5] (dot-dashed lines). In 
black we show the entire TNO sample. The most prob- 
able mass in TNOs for the combination of all surveys 
is Mtno = 0.020lQQQ3Af® while for our survey alone 

we obtain Mtno = 0-025to;oo7^e- These are consis- 
tent with each other a nd with the previous estimate by 
IBernstein et al.l 1)20041 ). This is not surprising since most 
of the mass is present in TNOs with sizes comparable to 
the size at which the distribution breaks. The slight over- 
abundace of TNOs in our survey with respect to other 



surveys yields a higher mass for the TNO population. It 
is important to note that in equation [Tl] the total mass 
diverges if either ai < 0.6 or 02 > 0.6. We also see in 
Figure [16] that for the results of our survey alone there is 
a long tail to higher masses. This is due to the poor con- 
straint on the bright end of the TNO luminosity function 
given the limited areal coverage of our survey (2.83 deg^) . 
However, the combination of all surveys yields a better 
constraint, and we obtain convergent masses for all steps 
in our MCMC run. 

When we consider the Classical and Excited popula- 
tions separately the mass distributions change. In Fig- 
ure [in] we show the mass in Classical objects in green 
and that in Excited objects in red. Using all the surveys 
the mass in classical objects is very well constrained to 
be Mela = 0.008 ± O.OOlAf^. Based on our survey alone, 
we find Mela = 0.013 ± 0.003Afe. The overabundance 
of Classical objects in our survey is responsible for that 
seen in the entire TNO population. 

The mass in Excited objects using all surveys is M^xc = 
O.O1O1q'qq3M0, larger than that found for the Classical 
TNOs and is also less well constrained, with a long tail 
to higher masses. This reflects the relatively poor con- 
straint on the size distribution of Excited objects, where 
the limits are set by what values for the exponent of the 
power law size distribution are considered to be physi- 
cally plausible. With only 18 Excited objects in our sur- 
vey we have a very poor constraint on the individual DPL 
parameters. However, the mass is well constrained. We 
find Mexc = O.OO5tQOQ3M0, less than the mass in Clas- 
sical TNOs. This is due to the relative under-abundance 
of Excited objects in our survey. This can be explained 
by the fact the survey was conducted in the direction of 
Uranus, separated about 18.5 deg from Neptune, where 
we expect Plutinos to be near apocenter and hence faint 
and under-represented. 

6. CONCLUSIONS 

We have presented a TNO survey that is both deep 
(i?5o=25.6) and broad (-2.8 deg^), finding 82 TNOs. 
The survey is very well characterized and simple, reach- 
ing its limiting magnitude in single exposures. 

We have studied the luminosity function of the TNOs 
in our survey. We found a significant deviation from a 
single power law behavior in the cumulative function at 
i?~25. We have shown that our data are consistent with 
a single power law, and with many other shallower sur- 
veys, if we consider only objects brighter than i?=24.5. 
We have also demonstrated that the apparent deviation 
from a single power law is not an artifact of our detection 
efficiency. 

Whether our data support a break in the luminosity 
function is a matter of statistical analysis. We compared 
two models, one where the distribution increases expo- 
nentially with a single power law and one where there 
are two different slopes in the sampled magnitude re- 
gion, and compute the t otal pro bability of each model 
with Bayesian statistics ([Gregory ,2005.) (See details in 
Section 15. ip . The ratio of the total likelihood for a dou- 
ble power law and a single power law model is ~26. This 
can be interpreted as the DPL model being 26 times more 
probable than the SPL given our data set. 

We conclude that our survey provides significant evi- 
dence for a break in the TNO luminosity function. This 
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is the first survey that is able to make such a claim with- 
out relying upon the results of other surveys. Our result 
is easy to interpret since we do not have to make as- 
sumptions about the distribution of objects in different 
parts of the sky. Nonetheless, the comparison with other 
surveys is fundamental since there are published searches 
that sample the same magnitude region. We have consid- 
ered most of the pubHshed data up to July 2007 regarding 
su rveys of the trans-nep tunian space in the same spirit 
of iBernstein et all |2G04). Again, our double power law 
model accurately describes the cumulative number den- 
sity for all surveys combined. 

Only two ground based surveys are as deep as the 
present, and they have not seen a signi ficant deviation 
from a single power law. The survey of iGladman et all 
(|2001;) covered much less area and, consequently, discov- 
ered many fewer TNOs in this magnitude range (17 ob- 
jects for the entire survey). Given the sm all numbers, our 
result s are not incon s istent with those of IGladman et all 
(|2001h . iFraser et all ()2008D report the combined results 
of surveys taken at different ecliptic latitudes and longi- 
tudes. They fit for a single power law but account for 
variations in the sky surface density, that may be due to 
surveying at different ecliptic longitudes and latitudes, by 
allowing an offset in the luminosity function zeropoint for 
each survey. This substantially increases the number of 
free parameters and, we believe, allows deviations from 
a single power law within individual surveys to be ob- 
scured when the results of several surveys are combined. 
We believe that this explains th e difference bet ween the 
present results and those of Era ser et al.l (j2008f ). 

We make the assumption that all objects are located 
at the same distance, so the luminosity function can be 
translated into a size distribution. For every object with 
a reliable distance estimate a nominal size can be com- 
puted (we assume an albedo of p—OM). The size dis- 
tribution of our survey was compared with the single 
distance approximation and we showed they agree. We 
then interpret the DPL size distribution. 

The break in the size distribution reflects the size at 
which coUisional processes take over gravitational ones. 
This is, the largest object that is expected to be dis- 
rupted in a collision in the age of the solar system. 
The best DPL model for our survey features a break 
at D=130 (p/0.04)-°-5 km bodies while for all sur- 
veys it is at D=100 (p/0.04)~°-^ km. Current mod- 
els expect the break t o occu r at smaller sizes, D < 
50 km for iPan fc Saril (l200l and D < 100 km for 
iKenvon fc BromlevI (|2004f K We consider these models 
to be consistent with our result given the assumptions 
on poorly constrained quantities like the albedos on the 
observational side as well as initial conditions in the the- 
ory are not well constrained. The effect of a distribution 
of albedos and a possible correlation with object size and 
heliocentric distances should be studied. 

The inclination distribution for our survey is consis- 
tent w ith what is expected from previous results (Brown 
l2001h . However, we do not have enough objects to do 
a detailed study of the distribution. We do, however, 
separate our population in classical ( "cold" ) and excited 
("hot") objects. We study the size distribution of these 
samples a nd find them to s how d ifferences as done pre- 
viously bv IBernstein et al.l (|2004l ). 



We calculate the probability distribution for the to- 
tal mass in TNOs, Classical and Excited objects that 
are consistent with our observations and all considered 
surveys. For all surveys combined we find Mtno = 
O.O2Olo;oo3M0. It is interesting to note that for the 
classical population the mass is very well constrained 
to be Mela = 0.008 ± O.OOIM® while the excited pop- 
ulation gives a larger and poorly constrained mass of 
Mexc = O.OIOIq'oos-Mq. This provides evidence for a dif- 
ference between the "hot" and "cold" populations. Our 
survey gives a consistent but slightly higher answer for 
classical objects that we believe is due to the local over- 
abundance of objects in our survey. We only have 18 ex- 
cited objects in our sample, too few to constrain the pa- 
rameters of the luminosity function, but enough to show 
there is an under-abundance of excited objects in our 
survey. This is explained by the direction of our fields, 
close to where most of the Plutinos come to apocenter. 

Given th e size distribution we ca lculate a distance bias 
correction (Trujill o fc Brownll2001h . We then obtain the 
real distance distribution of objects, assuming we are 
just as likely to find faint objects that are close as those 
that are far. Our survey is very well suited to detect- 
ing objects that show slow parallactic movement (dis- 
tant); our detection efficiency is essentially independent 
of rate for rates larger than 0.9 'Vhr (di s tances closer 
than 1 50 A U). According to iDqnesi ([1997| ). iJewitt et all 
(|1998[ ). and iTruiillo fc Brownl (|200lD ~the fraction h of 
objects found outside 48 AU should be about 40% for 
a population with a smooth brigthness distribution that 
extends beyond 50 AU. In our sample there are 73 TNOs 
with reliable distance estimates, of which 71 are located 
between 30 AU and 47 AU, and only two at ~50 AU, 
accounting for h=3%. Once we take into account the bi- 
ases associated with distance these numbers indicate an 
abrupt drop in the radial density of the Kuiper Belt. If 
we also consider the size distribution break found in our 
sample we also rule out the existence of a far population 
of TNOs near the plane of the ecliptic. We have found 
more evidence for an edge of the Classical belt population 
at around 47 AU and placed a constraint on the surface 
density of objects for an unseen population at 60 AU of 
8% that of the observed Classical Belt. We also set a 
minimum distance for a "belt-like" population with the 
same mass as that of the Classical Belt of 110 AU. 

Deeper surveys will help better constrain where the 
break in the luminosity function occurs and complete 
the picture of the trans-neptunian space. The size distri- 
bution would be better determined if these surveys are 
also careful in obtaining foUowup observations to mea- 
sure accurate distances for faint objects. 
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TABLE 1 — Continued 



Name MJD RA Doc Rmag ^Rmag dRA/dt dDec/dt dpar dtari i Aa{Uranus) 
["/hr] ["/hr] [AU] [AU] [deg] ['] 

All 82 trans-ncptunian objects found. The second night data were used when possible. The measured magnitude in the R filter with nominal errors is shown in Rmag- ^Rmag is a model for 
the photometric error based on the measure magnitudes of inserted, synthetic objects. dYlA/dt and dDcc/dt arc estimates of the measured motion of the object. The distance dpar is calculated 

with the assumption of a circular orbit, dbari is the barycentric distance estimate and i is the inclination estimate given by the Orhfit code IjBernstein Sz Khushalanil|2000f) . /S.a{Uranus) is the 

a b 

projected distance to Uranus during the observations. These objects were not found in the second night of observations. The result is unconstrained. 
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TABLE 2 
Surveys* 



Paper Q R50 A^C^ A^E^ N^^s N^^p 

deg2 



Chiang fc Brown (1999) 
Gladman et al. (2001') '^ 
Truiillo e t al. (2001)'' 
Allen et a l. (2002) 
Bernstein et al. (2004) 
Petit et al . (2006) N 
Petit et a l. (2006) U 
Eraser et a l. (2008) 
This survey'l 



0.01 


27.0 


1 


1 


2 


1 


0.322 


25.9 


7 


8 


15 


15 


28.3 


23.7 


38 


27 


71 


64 


2.30 


25.1 


15 


15 


30 


39 


0.019 


28.5 


3 





3 


5 


5.88 


24.2 


6 


21 


27 


22 


5.97 


24.6 


16 


20 


36 


34 


3.0 


20.8 


36 


31 


67 


74 


2.83 


25.69 


54 


18 


82 


74 



Details of the surveys considered in this work. O is the to- 
tal surveyed area. Rqq defines the R magnitude at which the 
survey's detection efficiency is 50% its maximum efficiency. The 
total number of objects discovered that had magnitude brighter 
than that at which the survey is 15% its maximum efficiency 
is Nobs-, as defined in iBernstein et al.| | |2004D . The expected 
number of objects for each survey given our most likely DPL 
luminosity function model for all surveys combined (see Fig. 
|13l l is Nexp- ^ Objects with inclination z < 5 deg and at a 
distance 38 AU < d < 55 AU arc considered as Classical Nc 
and the rest as Excited Ne- Based on Table 3 and comments 
in lGladman et al.| ( |200lD .^ Based on Table 2 and comments in 
IBernstein et all 1120041) . In the Classical and Extended classifi- 
cation We only considered objects for which there was distance 
and inclination information. 



TABLE 3 

DPL Parameter Estimation* 



Survey 




"1 


02 


0'23 


Req 


All surveys 
This surveyj 


TNO 

Classical 

Excited 

TNO 

Classical 


n 7c;+0.l2 

^■'0-0.08 

1 4+0.1 

l-*-0.3 

n fii+o-07 

7+0-2 

1 9+0.3 
-'-•^-0.4 


-0.3«-i 


i.5o«:l« 

0.82t0-i3 

o.68«:g| 

n O+0.5 

^•'i-o.s 

1 c+O.S 
l-t'-0.5 


24.8i;5-5 

23.3tli 
25.7+°i 
24.31°:? 
23.6to,r 



Best fit parameters and 1-cr confidence limits based on MCMC simulations. 
All surveys arc detailed in Tablc[2] ^ In this survey there were only 18 excited 
objects, too few to constrain a 4-parameter model. However we could fit a 
SPL with Q^O.62 ± 0.12 and Ho^24.2 ± 0.3 to this population. 
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Fig. 1. — Detection efficiency as a function of magnitude, with an error given by the number of objects implanted and found in each bin. 
The fitted curve corresponds to Eq.[2l where the best fit values are A=0.88±0.01, iJ5o=25.69 ± 0.01, «)i=0.28 ±0.04 and «;2=0.88 ± 0.15. 
-R50 corresponds to the magnitude at which our method is 50% as efficient as its maximum detection efficiency. 
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Fig. 2. — Histogram of the fraction of objects recovered as a function of rate. Bins are chosen to have similar numbers of objects. This 
demonstrates that our detection efficiency does not depend significantly upon the rate of motion. 
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Fig. 3. — Histogram of the magnitude error (AR) as a function of R magnitude for all implanted objects. Tlie error is defined as the 
difference between the implanted and measured magnitudes for the synthetic population. The dashed line is a gaussian of width ~0.1 mag. 
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Fig. 4. — The error in magnitude for synthetic objects as a function of magnitude is shown for different magnitude bins. The error 
is defined as the FWHM of the best-fit gaussian to the histogram of errors for all objects in each bin. The error bars correspond to the 
calculated uncertainty of the FWHM. Tlie curve is a quadratic fit to the data and defines the error estimate used for AKmag in Table [T] 
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Fig. 5. — Rate of motion in the sky for every TNO. Objects observed only on one night only are represented by triangles. The ecliptic 
motion is overplotted as a solid line. 
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Fig. 6. — Contours of the SPL likelihood function. The maximum lik e lihood point is marked with a dot (o=0.51, iJo =22.6). Marked 
with a triangle is the best value f o r the parameters based on lPetit et al.l II2006I) . (a=0.76, -Ro=23.3). This shows the discrepancy between 
our result and that of lPetit et all pOOSl V 
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Fig. 7. — Contours of the SPL likelihood function for our sample limited to R SJ24.5. The m aximum l ikelihood point is marked with a 
dot (o=0.69, -Ro=23.0). Marked with a triangle is the best value for the parameters based onPetit et ak (20061). Both results consistent 
with each other. This shows that our survey agrees with previous surveys if we consider the only the range of magnitudes to which those 
surveys are sensitive. 
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Fig. 8. — Contours of the SPL likelihood function for our sample limited to R ^25.2, wher e our survey is 70% efficient. The maximum 
likelihood point is marked with a dot (a=0.57, -Ro=22.8) and the triangle is the IPetit et al.l l[2006) result. We see that both results are 
inconsistent at more than a 2-u level. This demonstrates that our result does not rely on the detection of objects at magnitudes where our 
detection efficiency is declining. 
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Fig. 9. — The Probability Density Function for a and Rq from the MCMC simulation is shown as a histogram. The Likelihood Function 
in Figure |6] shown as the marginal probability over each parameter is plot as the solid curve. The solid, heavy line indicates the global 
maximum obtained by the MCMC run and the thin lines indicate the l-cr credible region of the parameter, inside of which we find 68.3% 
of the probability. 
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Fig. 10. — The DPL likelihood for our survey as a function of all parameters is shown in each window. The most likely parameters and 
their 1-a confidence regions, represented by the solid, heavy line and the two thin lines, are oi=0.7_qj^, a2=0.3_Q'2i <'"23=2.0_q'5 and 



A Subaru Archival Search for Faint TNOs 

Size[km] 
^00 



23 



100.0 



V 



10.0 - 




Fig. 11. — The cumulative number density for our survey. The best previous model is plotted in the short-dashed line. Our most likely 
solution for the single power law is plotted in the long-dashed line. The best DPL fit is shown as a solid line. The quoted size corresponds 
to an object at 42 AU and 4% albedo. 
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Fig. 12. — The cumulative number density for all surveys in Table[2] The best previous mod el is plotted in the black dashed line. Our 
most likely double power law is plotted in the long-dashed line. The most likely DPL (see Fig. I13II co nsidering all surveys is plotted as a 
full line. The apparent bump in density at around _R~25.8 corresponds to 5 objects in lGladman et aT.i (i2001i ). 
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Fig. 13. — The DPL likelihood function marginalized over each parameter for all the surveys in Table |2] We see the maximum and 68% 
confidence region. The most likely value for each parameter and 1-cr confidence limits are: ai=0.75_Qgg, a2=0.23_g5^, 0-23=1. 50_q]^2 
and Req=2A.8+°ol. 
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Fig. 14. — Number of TNOs observed in 1 deg^ as a function of size. The solid line shows the model based on our survey. The dashed 
one is the model that considers all surveys. Both models are properly sealed to match the density observed in our survey. 
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Fig. 15. — Magnitude and distance for all 82 TNO's found. The black dots are objects observed in both nights and the triangles are 
those with only one night's observation. We assume a 4% albedo to plot the constant size curves for 100, 200 and 400 km in black. 
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Fig. 16. — Mass distribution for our population of TNOs. We are extending our result over a solid angle of 360 X 6 deg , assuming the 
fraction of the population that is within this solid angle / is 1. The solid lines represent the results from our survey alone while all surveys 
in Table [5] are shown as dot-dashed lines. The black lines correspond to the whole TNO sample; green and red are used for the Classical 
and Excited sub-samples, respectively. 
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Fig. 17. — Shape of the debiascd distance distribution of TNOs assuming a constant albedo=4%. The triangles assume a size distribution 
with a power law of index q=A. The dark points assume the DPL size distributions. The bias corrections /3 are overplot as a dot-dashed 
line for a single power law with exponent g = 4, a solid line for a broken power law based on our survey and as a dashed line for parameters 
based on all surveys combined. 
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Fig. 18. — Inclination probability distribution of TNOs in our survey. 



